Quantum deflation of classical extended objects 
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It is shown, that extended particle-hke objects should infinitely long collapse 
into some discontinuous configurations of the same topology, but vanishing mass. 
■ Analytic results concerning the general properties and asymptotic rates of such a 

Q^ \ process are given for 1+1-dimensional soliton models. 



O 



(N 



O 



hep-th/9410160 



At present, classical particle- like solutions (henceforth solitons) play a sig- 
nificant role in the QFT, since they are considered as reliable models of quantum 
> '. extended objects, e.g. existing baryons [1-3]. It is widely beheved, that if the 

classical soliton is stable due to either topological or dynamical reasons, then 
upon quantization by existing techniques its quantum descendant should most 
^ , likely reveal the same stability properties [1-3,6-9,11-15]. The purpose of this 

■ paper is to present arguments, that in the QFT, under rather general circum- 

stances, such a correspondence cannot take place. Namely, we'll describe a 
p peculiar quantum effect, that makes possible a special kind of collapse of such 
; extended objects into certain almost discontinuous configurations of the same 

topology, but vanishing mass. 
^ ■ Such configurations appear quite naturally within the explicitly lorentz-co- 

; variant treatment of the soliton Heisenberg field by means of Newton- Wigner 

position variables [4]. The latter emerge as the direct result of summing up 
the pertinent terms in the perturbative solution of eqs. of motion starting from 
the static classical solution [5] and serve as the covariant collective coordinates 
of the soliton, which restore the broken Lorentz symmetry and provide the 
cancellation of corresponding zero modes [6-9]. These covariant coordinates 
give rise to an effective dynamics, where all the differential operators of the 
initial theory are replaced by finite differences with the step, proportional to the 
(effective) Planck constant of the theory Tl [4,10]. In turn, such effective theory 
yields a series of c-number copies of classical solution, labeled by an integer 
A^, which are quite analogous to and reveal the same topological properties 
as the initial soliton, whereas their masses are much less than the classical 
one and vanish for h/N [4]. Within such a picture the soliton collapse 
appears quite naturally as a transition from the classical soliton state into its 
Nth copy with the same topological charge. The pecufiarity of this process 
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is, that the dynamics of the cohapse can be successfully understood in the 
quasiclassical approximation, leading to a universal power law for the energy 
loss into radiation in the asymptotics. Quantum origin of the effect manifests 
here rather in vanishing energy of the final soliton state and in the model- 
dependent tunneling effects, which determine the collapse probability Pcoii — 
ex.p{—S/h), where the value of S can be zero as well. 

To avoid unnecessary complications we'll consider these effects within the 
generic l+l-dimensional model of a nonlinear scalar field (p{x,t), described by 
the Lagrangean density 

Ci^) = '^id,^f-Ui^), (1) 

which possesses a classical solitary wave solution 

, ^ / X — vt \ , , 

^,ix,t) = u{-^=^j. (2) 

Soliton quantization with special emphasis on relativistic properties is imple- 
mented by means of the exact covariant center-of-mass coordinate [4-9] , which 
in l+l-dimensions is given by 

</ = lO + 4). (3) 

where L is the generator of Lorentz transformations and H is the total Hamil- 
tonian of the system. The operators H, L and the generator of the spatial 
translations P obey the Lorentz algebra 

[L, H] = ihP, [L, P] = ihH, [H, P] = 0. (4) 

The general analysis says, that the c.m.s. coordinate enters the Heisenberg 
operators always in the combination x — q{t), while the spatial translation is 
induced by q ^ q + a [5,7-9]. For our purposes it would be more convenient to 
introduce g in a dimensionless form by means of the substitution 

x^a^,t), (5) 

where 

^(x, t) = H^/\x - q{t))H^I'^ = Hx-Pt-L. (6) 

So in the vicinity of the soliton solution the Heisenberg field (p{x,t) is repre- 
sented as 

^ix,t) = fiax^t))+^a^,t),t). (7) 

In eq.(7) we take care of that the operator- valued argument of the soliton field 
might affect the soliton shape to be different from the classical one and therefore 



use new notation f{^) for it instead of that in (2). The meson field ^{^,t), 
considered as a function of c-number arguments, is taken independent of the 
generator P of spatial translations of the algebra (4), while the commutation 
relations between H, L and ^{^,t) are determined in such a way, that enables to 
provide the covariance of the whole field. The general procedure of such type has 
been considered in detail in refs. [6-9] . However, within the present approach 
we'll deal mainly with the soliton field f{^{x,t)), which is a covariant operator 
by construction. Note also, that ^{x, t) commutes with the mass operator = 
^2 _ p2^ gQ considered as a c-number by dealing with operators 

of the type f{^{x,t)). 

It can be argued, that the substitution (5-7) represents the general form 
of the Heisenberg field in the one-soliton sector. Indeed, it was shown in [5], 
that the replacement (5) appears as a direct result of summing up the perti- 
nent terms in the perturbative solution of equations of motion starting from the 
single classical soliton. It can be also shown [7-9], that within the expansion 
in inverse powers of the soliton mass this method is equivalent to the canoni- 
cal quantization via collective coordinate [11-13]. A quite diflterent relativistic 
framework for soliton quantization based on the BRST approach has been pro- 
posed recently in refs. [14,15]. However, the latter method turns out to be 
not quite appropriate to study quantum effects considered below, that are es- 
sentially non-analytic in h and so lie beyond the quasiclassical perturbation 
expansion. 

The main consequence of the substitution (5) is that in the resulting ef- 
fective theory all the differential operators of the initial theory are replaced 
by finite differences of special form. Namely, the dynamical equation, that 
determines the invariant soliton shape /(^), reads [4] 

/(^ + ^n) + /(^ - ih) - 2/(0 + ^v'ifiO) = 0, (8) 

where it is convenient to extract the dimensional (mass) parameter m from the 
interaction term and to introduce the dimensionless soliton mass 

U{ip) = m^V{(f), ii = M/m. (9) 

Then the eq.(8) is written completely in terms of dimensionless variables. (Re- 
call, that by means of a suitable redefinition of variables [1,4], in l-|-l-dimensions 
we can always treat h as the dimensionless expansion parameter for the soliton 
sector.) 

The mass of the soliton is determined from the following equation 

i'' = i^^h [^(^ + - -^(^^II/K - «) - /(«)] + jdi vmo) , (10) 



where the status of ± ih) is the same as in eq.(8). The detailed discussion 
of eqs.(8) and (10) is given in ref. [16]. 

The main consequence of the finite-differences in the dynamical equations 
is the appearance of what should be called quantum copies of classical solu- 
tions [4]. Namely, seeking the solution of eq.(8) in the form of the exponential 
(Dirichlet) expansion, we obtain the following series of approximate solutions, 
valid for h sufficiently small, namely 

fNiO = u{aNO, (11) 

where is an integer, 

1 27riV 

Q^TV = — + (12 
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A^iv = A^[/iv] is the corresponding dimensionless mass of the field, and a is equal 
to the mass of the elementary excitation of the field (meson) divided by m. 

For = the singular part m. disappears and the corresponding 
solution /o(0 ~ ""(^/a^o) survives the limit ^ 0. So /o(0 should be considered 
indeed as the quantum descendant of the classical solution, and so in what 
follows the case = will be referred to as the quasiclassical one. The masses 
of soliton copies are determined from eq.(lO) and for small h and 7^ are 
given by 
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\\aN\ ^|A^|V27r; 
while fiQ is the classical soliton mass 



(13) 



lj,Q = jdx u^{x). (14) 

The most essential point here is, that /ijv turn out to be 0{Vh) in magnitude, 
hence sufficiently smaller, than the classical mass /iq, but larger, than the typical 
meson energy, that is of order 0{h). Another peculiarity in the result (13) is 
that fiM decrease for increasing |A^|. These properties will be crucial for the 
soliton collapse. 

Now let us consider these and other properties of quantum copies using 
as an illustrative example the spontaneously broken (/?^-model, that is given 
by the selfinteraction potential V{ip) = | (1 — 9:?^)^, and yields the classical 
kink solution ipc{x) = tanhx. The nonclassical kinks are given by the shape 
functions 

/7v(0 = tanhttTv-C, (15) 

where = 1/hn + T^N/h. (Recall, that in the (/9'^-model with such potential 
the mass of the elementary meson is 2m, hence (7 = 2.) 



An important property of quantum copies is that for all 7^ they carry 
the same topological charge Q as the classical soliton. For h/N /n^O 
converge to the limiting shape fooiO the form of the discontinuous step- like 
classical solution of the same topology. In the case of an odd classical solution 
of the (f^-tj-pe, what includes the most important soliton models, we can write 
without loss of generality 

fNiO - fooiO = signe . (16) 

According to eq.(13), the mass of foo{0 vanishes. This result is essentially non- 
trivial, since in the purely classical case the kinetic energy of the discontinuous 
solution like signx is proportional to Jdx S'^{x) and diverges. The reason is, that 
the energy of quantum copies is determined by the operator fields /Ar(^(x,t)), 
rather than by the c- number shape functions /iv(^). In particular, for suffi- 
ciently small h (large N) the nonclassical kinks (15) can be approximated as 

MO ^ tanh (-^^j , (17) 

hence are i^-periodic, and so by substituting them into the differences [/(^ ± 
ih) — f{0] of eq.(lO) the latter vanish identically. Thus the quantum origin of 
soliton copies is crucial for such specific behaviour of their masses. 

On account of these properties of quantum copies we should expect, that 
any given soliton state will be unstable under transition with AN > 0, while 
emitting a corresponding number of mesons. Since the topological charge Q 
remains unchanged, it should be just the soliton collapse, rather than the soliton 
decay, that proceeds with AQ ^ 0. Now let us present the general analysis 
of transition from the quasiclassical soliton state into its Nth quantum copy. 
It would be natural to describe such process in terms of suitable dynamical 
variables, that continuously transform the initial configuration /o(C) = w(«oO 
into the final one /iv(0 = w(aivO- The simplest path without infinite potential 
barriers, that enables such transition, is given by the (linear) superposition of 
initial and finite soliton states, both taken in the rest frame, 

f{^.t)=Co{t)fo{0+CN{t)fNiO , (18) 

where the finite energy of the intermediate configurations is provided by the 
subsidiary condition, imposed on scaling parameters, namely 

co{t) + CNit) = 1. (19) 

So the path (18) reduces to a one-parameter family of trajectories 



fii, t) = /o(f ) + c(t)A/(f ) , 



(20) 



where 

A/(0 = fNiO - /o(0 (21) 
is a regular function with (exponentially) decreasing asymptotics. The collaps- 
ing process is given by the continuous evolution of c{t) subject to boundary 
conditions 

c(0) = , c{t) = ; c(oo) = 1 , c(oo) = 0. (22) 

It would be natural to assume, that the collapse is an adiabatic process 
(what is confirmed by explicit calculation below). Then the energy loss into 
radiation will be very low, hence $(^,t) will be small compared to f{^,t). 
Granted this, the meson field can be treated perturbatively in the slowly varying 
c-number soliton background f{i,t), while the back reaction of ^{^,t) on the 
dynamics of c(t) being neglected. So in the first step we have to consider the 
proper dynamics of c{t) without radiation. For these purposes we insert the 
collapsing soliton field (20) into the initial field Hamiltonian and evaluate the 
energy (mass) of the configuration as an explicit function of c{t) and c{t). The 
result looks like eq.(lO) by adding the terms with c{t), namely 

l,'='^c\t)ld^ [A/(Of + ^|{c(^)/ciC [il-cosnd^)fi^,t)] A/(0} + 

+ jdin[i,t), (23) 

where 

n{i. t) = ^ { [(1 - cos nd^) /(^, t)f + [(sin nd^) /(^, t)f} + v(/(c, t)). (24) 

Note, that by integrating over in eqs.(23) and (24) we can always take 
the final configuration as the discontinuous /oo(0- definitely correct for 

transition to A/" — ^ oo, but for h small such replacement would be valid actually 
for any transition from the quasiclassical state into N 0. 

Calculation of the finite-difference kinetic term T in Jd^ ^{^,t) proceeds 
by means of the following relations [4] , valid for any Nth soliton copy 

isinhd^) fNiO - —u'iaNO. (1 - coshd^) fNiO - \ ^u'\aNO , (25) 

where the prime denotes the derivative with respect to the argument, and the 
residual in both expressions (25) is 0{fi). Therefore, for h small the term 
[(1 — coshd() /n] can be certainly neglected while compared to [(sin/ic}^) f^] . 
Then on account of the eq.(lO) we obtain for the kinetic term of the eq.(24) the 
following expression 

[ismnd^)fi^,t)f = n\cl + c%) + 2 coCN Id^u'iaoOu'iaNO- (26) 



The last term in eq. (26) turns out to be 0(/i^/^), hence in \/h times less than 
the first one, and so can be dropped to the lowest order. So the answer for T 
acquires a very simple form 
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r = ^(c^ + 4) = |[c^ + (l-c)^]. (27) 

Applying the same procedure to the second term in the eq.(23), we can estimate 
it as 0(/i^), and so for an adiabatic process with bounded derivatives of c{i) 
this term can be neglected as well. So finally the resulting expression for the 
energy of the collapsing soliton configuration (20) can be written in terms of 

c(t) and c(t) as 

.2 _ 7^2 , 

2 

where 

l = \di [A/(0]' , W{c) =\diV [/o(0 + cA/(0] • (29) 
For the (^^- model the straightforward evaluation of integrals in eqs.(29) gives 



,2 



= T^^' + ^ V + (1 - + ^(^) ' (28) 



7 = ?(21n2-l), W[c) = ^-{c-\f 



16\ , 2 2 
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(30) 



Now let us consider more carefully the most general features of the be- 
haviour of c(t), that are irrespective of the specific properties of the model. 
Firstly, by definition W{c = 1) = 0, hence the point c = 1 corresponds to the 
absolute minimum of W{c), while W{c = 0) = So the point c = 1 wih 

always lye within the classically allowed region for the motion of c(t), whose 
velocity at this point is = fi^/'y. 

Further, with respect to the shape of W{c) on the interval [0, 1] one has to 
distinguish between two different situations. The first one concerns the case, 
when both the points c = and c = 1 belong to the same classically permissible 
region of motion. Then c{t) will swing nonlinearly between the neighboring 
turning points, that are given by Ci = and C2 > 1, while the point c = 1 will 
be the intermediate point of the trajectory. Indeed such a situation takes place 
in the (/?^-model. 

The second case corresponds to the situation, when W{c) yields a finite 
potential barrier between c = and c = 1 . Now in order to reach the final point 
c = 1 while starting from c = 0, we have firstly to penetrate the barrier, what 
gives the quantum-mechanical tunneling amplitude Ptunn for the probability of 
the whole collapsing process. As soon as we penetrate the barrier and so reach 
the classically permissible region, which the point c = 1 belongs to, there starts 
the process of nonlinear oscillations between the points < ci < 1 < C2, similar 
to described above. Since the tunneling proceeds at the energy close to the 



initial classical soliton mass /ig, that is large enough to use the quasiclassical 
picture, the tunneling amplitude can be estimated through the conventional 
WKB-approximation 

Ptunn — 6 ^ f (3-'-) 

where S is the corresponding classical action for the forbidden region. 

The most strict predictions can be made for the asymptotic behaviour of 
the collapse, when n is sufficiently small already and so the oscillations shrink 
to a small neighborhood of the minimum of VK(c), i.e. to the point c = 1. Then 
we can replace the second term in the eq.(28) by //^/2, simplifying it up to 

y = + W(c). (32) 

Now let us write 

c{t) = 1 + c{t) , (33) 
where c{t) is small, and represent the collapsing configuration as 

/(C,t) = /oo(0+c~(t)A/(C). (34) 

Since T^(/oo) = ^'(/oo) = 0? from the eq.(34) we observe at once, that the 
expansion of W{c) in powers of c is equivalent to the expansion around the true 
vacuum value, what gives 

(35) 

As a result, the limiting oscillations will be harmonic with the frequency 

v = G . (36) 

Such peculiar behaviour of the asymptotical dynamics of the system gives 
rise to one more general feature of the collapsing process. Let us assume, that 
there exists a set of discrete modes (bound states) 0n(a^) in the spectrum of 
meson excitations around the soliton. Then it might seem, that the transition 
from the initial to finite soliton states might occur through the excitation of 
certain However, actually such a picture doesn't take place, since the eigen- 
frequencies of the discrete spectrum lye always lower, than the meson mass, 
namely < < cr, and so due to the relation (36) the resonance in the me- 
son excitation will be just at the threshold of the continuum. Therefore, the 
excitation of the bound states will be suppressed compared to the creation of 
wavepackets, which remove the energy from the soliton. Moreover, it can be 
argued [16], that the mesonic modes around the final soliton state fall out of 
the dynamics at all. 



7 



To estimate the energy loss into radiation we have to look at the dynamics 
of the meson component t) during the collapse. We start with the obvious 
statement, that the energy of the full Heisenberg field (p{x,t), that includes 
both the collapsing soliton and the meson field, certainly conserves. It means, 
that if the starting point of the collapse is c = 0, c = 0, 11 = 110, the total mass 
of the field remains equal to /iq throughout the process. In turn, in the rest 
frame the wave operator Oxt acquires in terms of finite differences (after suitable 
redefinition of the time t mt) the following form 

Uxt^d^t^2^{cosnd^-l). (37) 

As a result, the initial field equation, whence the dynamics of $(^, t) should 
have been determined, can be written as 



d'i + 2^{co^hd^-l) 



^{e„t) + V'm,t)) = Q . (38) 



In the eq.(38) <p{^,t) stands for the Heisenberg field transformed to the ^- 
reference frame 

ifix^t) = 0i^ix, t),t), (39) 

whereas during the collapse 

^{U) = f{^.t) + ^U)- (40) 

The direct consequence of such initial conditions is that /o(0 satisfies the 
eq.(38) with the required precision. Therefore the region |c| <C 1, where f{^,t) 
is close to /o(0) cannot significantly contribute to the radiation per period 
of oscillations. The dominant contribution originates from the neighborhood 
of c = 1, where f{^,t) ~ fNi^,t)- During the initial stage of the process 
such a picture will hold, because fniO doesn't satisfy the eq.(38), whereas 
for t large the oscillations of c{t) will shrink to a small neighborhood of c = 1. 
There remains, of course, an intermediate region, where the oscillations won't be 
harmonic yet, while their amplitude will be large enough for the nonlinear effects 
to be significant. However, in the asymptotics the radiation will be indeed very 
small, while at the beginning of the process the nonlinear terms cannot play 
any significant role, since it takes some definite time for any nonlinearity to 
develop [17]. Therefore these effects can be essential only for a finite interval 
of time, and so cannot seriously affect the main properties of the process, while 
depending to a high degree on the specific features of the model. So in order 
to retain the generality we'll omit this nonlinear region. 

Now let us define more correctly the status of the meson excitations during 
the collapse. First of all, since the meson field must compensate the c-number 



difference in soliton shapes Af{^), to the lowest order t) should be treated 
as a c-number coherent wave. Further, any mesonic excitation should be re- 
ferred to the initial soliton state with N = 0, what we explicitly take account 
of by putting 

$(^,t) = l>(ao^,t). (41) 

Making use of relations (25), we can now easily verify, that for such a function 
the finite-difference wave operator (37) can be replaced by 



^^,t)^nytHy,t), (42) 



where y = ao^. Actually, the last relation is nothing else, but the manifestation 
of the quasiclassical nature of the sector with A^^ = 0. 

After these preliminary remarks we insert the decomposition (40) into the 
eq.(38) and linearize it with respect to ^{^,t) in the vicinity of /7v(0' while 
supposing, that the main contribution to $(^, t) appears from the neighborhood 
of c = 1. Then we get the simplified equation 

[ Oyt + Vih) ] !> = i, (43) 

where the radiation generating current j is determined as the residual of the 
eq.(38), applied to the collapsing soliton field in the vicinity of fNiO 

-M,t) = c(t)A/(0 + 2^ [(cos;.^^ - 1) MO] + v'iMO]- m 

By means of the relation (25) the difference term in j{^^t) can be represented 
as 

2 ^ [{coshd^ - 1) MO] ^ -4 ^"(^nO- (45) 
It is easy to verify, that for /i/A ^ we can write 

4 «"(«ive) - ^ (46) 

where Q is the topological charge and so doesn't depend on ?i. The factor a'^ 
in the r.h.s. of the eq.(46) survives the integration over with the smooth 
function of the type (41), but vanishes for h/N ^ 0. So the difference term in 
the radiation current can be dropped. 

Finally, we replace fniO the eqs.(43) and (44) by the hmiting /oo(0' 
for which V'{foo) = 0, V"{foo) = cr^. So in the lowest approximation the eq. 
of motion for the meson field reads simply 

(pyt + a')^y,t) = -c{t)Af{0. (47) 



Now, to avoid the excess of precision we retain in A/(,^) only the first term in 
the asymptotics for |,^| — > cx). Then for the case of an odd classical soliton we 
can represent A/(^) explicitly as 

A/(0 - e-^'^'sign?/, y = a^t (48) 

Motivated by the same reasons, in the last step we keep only the main harmonics 
of c{t) with the frequency z/ = cr, approximating the oscillations by 

c{t) = 1 + A cos ut. (49) 

Then the eq.(43) can be easily solved, whence for the radiation energy lS.S per 
the main period of oscillations At = 27r/z^ we find 



^ ^ ^ (Av^) f dk . . [At r-^ ^ 

A£: = 4^ ^ / — — sin2 —VP + ^2 



= 3.926 . (50) 

TT 



Thus, the radiation energy per period can be estimated as AS = Const x 
A^, where A is the magnitude of oscillations around c = 1. From this relation 
the asymptotic rate of the collapse can be easily deduced while noticing, that 
for sufficiently large t the remaining soliton mass is determined through the 
eq.(35), what gives 

= 7 u'A^ (51) 

Taking the average over the period, we can write down now the following equa- 
tion for the energy loss 

^ = -AS = -Const X fi^ , (52) 

where 

Const = (53) 

Therefore the asymptotical law for the remaining soliton mass /j,{t) reads 

Const 

Kt) = , (54) 

while the amplitude of oscillations around c = 1 and the radiation power de- 
crease as 1/t and correspondingly. 

So the soliton collapse turns out to be indeed a very slow process with 
the time scale given by the constant (53). It should be stressed, however, that 
this constant plays here a quite different role compared to the decay time r 
of the quantum-mechanical quasistationary state, that decays according to the 
exponential law exp(— t/r). In the latter case r determines the time interval. 



during which practicahy ah the excitation energy should have been emitted, 
whereas the soHton collapse reveals a slowly decreasing tail of radiation. 

The asymptotic rate of the collapse doesn't depend explicitly on the differ- 
ence parameter h (or h/c, when the speed of light c is restored explicitly). The 
reason is, that the underlying quantum origin of the effect manifests here in the 
vanishing energy of the final soliton state, whatafter the whole process can be 
successfully formulated within the quasiclassical picture. However, there exists 
another possibility for h to enter the resulting dynamics, namely through the 
collapse probability Pcoii- When there are no barriers, as for the (/?^-theory, the 
whole collapse proceeds in the real time, thence 

Pcoll = 1. (55) 

On the contrary, in presence of a barrier one has firstly to penetrate it, what 
gives rise of the tunneling amplitude Ptunn-, and so the collapse probability can 
be estimated as 

Pcoll = Ptunn ^ e-^/^ , (56) 

although in both cases the asymptotical behaviour is governed by the same 
power law (54). So in presence of tunneling the collapse acquires an explicit ex- 
ponential dependence on h, what apparently demonstrates the non-perturbative 
origin of the effect. Remarkably enough, this time h comes merely from the 
purely quantum-mechanical considerations, hence it has nothing to do with 
finite differences and so with the speed of light. So it turns out, that while 
emerging as an essentially relativistic effect, the collapse actually survives the 
non-relativistic limit c — oo. 

To conclude let us note, that such a collapse turns out to be an immanent 
feature of the particle-like classical solutions in the relativistic QFT. The rea- 
son is, that the energy density of such extended objects is locahzed, therefore 
the covariant c.m.s. coordinate is a well-defined quantity (see [18-20] and refer- 
ences therein). Then the finite-difference structure of the effective dynamics in 
covariant coordinates implies, that upon quantization the mass of the limiting 
step-like configuration /oo(0 ^^^^ vanish, rather than diverge. In this way the 
extended particle-hke objects induce their collapse by themselves. 

Actually, the derived asymptotic law (54) should be valid in more spa- 
tial dimensions as well, since the latter maintain the difference structure of 
the effective dynamics in covariant coordinates [10]. Moreover, the structure 
of corresponding difference operators coincides with those of the quasipotential 
approach on the mass shell in momentum space [21], what once more empha- 
sizes the relativistic origin of the effect. So we should expect, that the main 
features of the effect will survive the transition to more spatial dimensions. The 



application of the general framework developed in [10] to this problem will be 
reported separately. 
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